Unveiling New Magnetic Phases of Undoped and Doped Manganites 



Takashi Hotta 1 , Mohammad Moraghebi 2 , Adrian Feiguin 3 , Adriana Moreo 2 , Seiji Yunoki 4 , and Elbio Dagotto 2 
1 Advanced Science Research Center, Japan Atomic Energy Research Institute, Tokai, Ibaraki 319-1195, Japan 
2 National High Magnetic Field Laboratory, Florida State University, Tallahassee, Florida 32306 
3 Department of Physics and Astronomy, University of California at Irvine, CA 92697 
4 International School for Advanced Studies (SISSA), via Beirut J h 34014 Trieste, Italy 

(Dated: February 1, 2008) 

Novel ground-state spin structures in undoped and lightly-doped manganites are here investigated 
based on the orbital-degenerate double-exchange model, by using mean-field and numerical tech- 
niques. In undoped manganites, a new antiferromagnetic (AFM) state, called the E-type phase, is 
found adjacent in parameter space to the A-type AFM phase. Its structure is in agreement with 
recent experimental results. This insulating E-AFM state is also competing with a ferromagnetic 
metallic phase as well, suggesting that large magneto-resistant effects could exist even in undoped 
Mn oxides. For doped layered manganites, the phase diagram includes another new AFM phase of 
the C x Ei_ x -type. Experimental signatures of the new phases are discussed. 
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In the recent decade, the study of manganites - materi- 
als that show a remarkable Colossal Magneto- Resistance 
(CMR) (J - has been one of the most important areas 
of research in condensed matter This CMR effect 
occurs when the manganite ground-state changes from 
insulating to ferromagnetic (FM) metallic after a small 
magnetic field is applied. Based on the concept of two- 
phase competition |2j, the CMR behavior has been suc- 
cessfully qualitatively reproduced in computational sim- 
ulations employing resistor- network models 3]. However, 
more work remains to be done to fully understand Mn- 
oxides, both regarding their unusual magneto-transport 
properties and the nature of the many competing phases. 

The appearance of the FM metallic phase in man- 
ganites is usually rationalized by the so-called double- 
exchange (DE) mechanism, based on a strong Hund cou- 
pling between mobile e g electrons and localized t2 g spins. 
On the other hand, the insulating phase in manganites 
occurs due to the coupling between degenerate e g elec- 
trons and Jahn- Teller (JT) distortions of the MnOg oc- 
tahcdra, leading to the various types of charge and/or 
orbital orders observed experimentally 0. 

The parent compound of CMR manganites is undoped 
RMn03, where R denotes rare earth ions. For R=La, 
as is well-known, the A-type antiferromagnetic (AFM) 
phase appears, with the C-type ordering of (3x 2 — r 2 )- 
and (3y 2 — r 2 )-orbitals Q|. By substituting La by alka- 
line earth ions such as Sr and Ca, holes are doped into 
the e g -electron band and due to the DE mechanism, the 
FM metallic phase appears, with its concomitant CMR 
effect. Most of the discussion in manganites has centered 
on the many phases induced by doping with holes the A- 
type AFM state, at different values of their bandwidths. 
In this framework, it is implicitly assumed that the un- 
doped material is always in an A-type state. However, 
quite recently, a new AFM phase has been reported as the 
ground-state in the undoped limit for R=Ho @, 0] . This 
phase is here called the "E-type" spin structure follow- 
ing the standard notation in this context. It is surprising 



that a new phase can be still found even in the undoped 
material, previously considered to be well understood. In 
addition, the nature of the states obtained by lightly dop- 
ing this E-phase is totally unknown, and new phenomena 
may be unveiled experimentally in the near future. Over- 
all, this opens an exciting new branch of investigations in 
manganites since novel phases appear to be hidden in the 
vast parameter space of these compounds. A clear exam- 
ple has been recently provided by the prediction of a FM 
charge-ordered (CO) phase at x=l/2 SB, which may 
have been found experimentally already |9| . 

In this Letter, based on the orbital-degenerate DE 
model coupled with JT distortions, the ground state 
properties of undoped manganites are analyzed by us- 
ing mean-field (MF) calculations and Monte-Carlo (MC) 
simulations. In our phase diagram at x=0, the E-AFM 
phase is found to exist in a wide region of parameter 
space, adjacent to the A-AFM phase in agreement with 
the experimental results. The E-AFM phase is robust 
when the dimensionality and/or electron-phonon cou- 
pling are modified, and its strength originates in the in- 
trinsic nature of orbital-degenerate DE systems. In the 
phase diagram and at small electron-phonon coupling the 
E-AFM insulating phase is located just next to a FM 
metallic phase [l(|, suggesting that CMR effects could 
be found even at x=0. Light-hole doping x of the E-type 
phase is also discussed and another novel magnetic phase, 
defined as the "C x Ei_ x " phase, is found. The ubiqui- 
tous phase-separation tendencies observed when insulat- 
ing and metallic phases compete is also expected near the 
FM-C x Ei_ x boundary. 

The Hamiltonian studied in this paper is 
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where di aa (dibo-) annihilates an e g -electron with spin 
a in the d x 2_ y 2 (d 3z 2_ r 2) orbital at site i, and a 
is the vector connecting nearest-neighbor (NN) sites. 
The hrst term is the NN hopping of e g electrons with 
amplitude , between 7- and 7'-orbitals along the 
a-direction: it 
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t| a =t| b =t^ a =0 for a=z. Hereafter, t is taken as the en- 
ergy unit. In the second term, the Hund coupling Jh(>0) 
links 



electrons with spin s ; = J2 ia p d\ 1 a (T ^ d ^0 
(<x=Pauli matrices) with the localized i2g spin Si as- 
sumed classical with |Si| = l. Jjj is here considered as 
infinite or very large. The third term is the AFM 
coupling Jaf between NN t2 E spins. The fourth term 
couples e g electrons and MnO@ octahedra distortions 
A is a dimensionless coupling constant, Qu is the 
breathing-mode distortion, Q21 and Qs\ are, respectively, 
2 )-type JT-mode distortions, p\= 
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Ysa(<haadia<j -dl ba d iba ). The fifth term is the usual 
quadratic potential for adiabatic distortions and /3 is the 
spring-constants ratio for breathing- and JT- modes. In 
actual manganites, /3«2. In undoped manganites, all 
oxygens are shared by adjacent MnOg octahedra and 
the distortions are not independent, suggesting that the 
cooperative effect is even more important than for the 
doped case x>0. To consider this cooperation, here oxy- 
gen ion displacements are directly optimized |l2(. 

Let us first describe our two-dimensional (2D) results, 
since the essential physics behind the stabilization of the 
E-type phase can be grasped by using MC simulations 
with relatively short CPU times. The phase diagram on 
a 4x4 lattice is shown in Fig. 1(a). The continuous curves 
are obtained by comparing the energies of the compet- 
ing phases in the MF calculations, while the circles are 
obtained by monitoring the nature of the dominant spin 
correlation 5(q) in MC simulations. A typical result for 
5(q) is shown in Fig. 1(b). A new regime characterized 
by q=(7r/2, 7r/2) is clearly observed between the FM and 
G-AFM phases. The good agreement between MF and 
MC results in the region of interest shows the high accu- 
racy of the present MF calculations for manganites 0] . 

In Fig. 1(c), the spin and orbital structure of the novel 
intermediate phase (E-phase) is shown. Along the zigzag 
chains, <2g spins order ferromagnetically, but they are 
antiparallel perpendicular to the zigzag direction. In the 
three-dimensional (3D) case, the MF study shows that 
the pattern Fig. 1(c) just stacks along the z-axis, while 
the spin directions are reversed from plane to plane. Note 
that the orbital structure is the same as that of the A- 
AFM phase, namely, the staggered pattern of (3a; 2 — r 2 )- 
and (3y 2 — r 2 )-like orbitals. Our investigations show that 
the E-phase is robust at weak- and intermediate-A, but 
for A>1.5, the E-type regime narrows. 

A surprising aspect of our results is that the E-type 
spin arrangement is the ground-state for a wide range 
of Jaf, even at A=0, indicating that the coupling with 
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FIG. 1: (a) Phase diagram obtained using the 4x4 lattice. 
Open circles denote MC results and solid lines are obtained 
with MF calculations, (b) Spin correlation S(q) vs. Jaf at 
A=1.5. (c) Schematic view of the spin and orbital structure 
of the E-AFM phase. Solid and open symbols denote orbitals 
with up- and down-spins, respectively. Note the presence of 
zigzag chains for each spin orientation, (d) Density of states 
(DOS) obtained by MC simulations at A=0.25. Solid and 
dashed curves denote the results for E- and FM-phases. (e) 
Gap at the Fermi level for the FM phase vs. A. The mag- 
nitude is evaluated from the DOS obtained by MC simula- 
tions, (f) Total energy vs. magnetic field for E-AFM and FM 
phases, (g) Phase diagram in the (Jaf, T) plane at x=0 and 
A=l, using a 4x4 lattice. Estimations of the Curie and Neel 
temperatures are evaluated using MC simulations (see text). 



JT phonons is not a necessary condition for its stabiliza- 
tion. This is in sharp contrast to the case of the A- AFM 
phase. To understand this point, it is instructive to study 
the e g -electronic structure of the zigzag FM chains that 
appear in the E-phase spin arrangement. Taking Jh as 
infinity for simplicity, the e g electrons move only along 
the zigzag FM chain, and cannot hop to the adjacent 
FM chains. The dispersion energy for e g electrons in this 
zigzag FM chain is given by £fc=(2/3)(cos fc±V cos 2 k + 3) 
and (2/3)(— cos fc±Vcos 2 k + 3), indicating that there 
appears a large band- gap equal to 4t/3 at half-filling. In 
fact, even on the 4x4 cluster, there is a clear gap of the 
order of t for the E-phase in the density of states (DOS) 
(Fig. 1(d)). Since t*=— for n^v, the sign in the hop- 
ping amplitude changes periodically in the zigzag situa- 
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FIG. 2: (a) Energies of the FM, A-AFM, E-AFM, C-AFM, 
and G-AFM phases on 4x4x2 bi-layer lattices (A=1.5). Solid 
circles indicate the results of optimizations, while lines denote 
the MF results, (b) Phase diagram for the 4x4x4 cubic lat- 
tice. All solid lines emerge from MF calculations. 



tion, leading to a periodic potential for e g electrons and 
its concomitant band-insulator nature. In other words, 
the E-type phase is stable due to the zigzag geometry of 
the FM chains that induce a band-insulator |15|. 

Another related interesting point is the orbital struc- 
ture of the FM phase. In the strong-coupling region, an 
orbitally ordered (OO) state appears, essentially with the 
same pattern as that observed in the E-AFM phase dis- 
cussed above. All MnOg octahedra are distorted at x=0 
and the cooperative effect is essential to determine the 
OO pattern for large A, irrespective of the spin struc- 
ture. However, when A decreases OO disappears and an 
orbital disordered (OD) new phase is observed. This is 
a metallic phase according to the DOS in Fig. 1(d). The 
OO-OD transition is monitored by the gap size at the 
Fermi level in the DOS (Fig. 1(e)). The transition ob- 
served using a 4x4 cluster is robust, although the actual 
critical value may change in larger systems |l6j . 

Note that the OD/FM phase is next to the insulating 
E-AFM phase for A<0.5 (the E-phase is insulating with a 
DOS gap both at small and large A). Such a result is new 
in the study of undoped manganites. Since the competi- 
tion between FM metallic and insulating phases is at the 
heart of the CMR phenomena, by tuning experimentally 
the lattice parameters in RMn03 it may be possible to 
observe the magnetic-field induced metal-insulator tran- 
sition even in undoped manganites. In fact, small mag- 
netic fields of a few Teslas indeed induce a phase tran- 
sition between E- and FM phases (see Fig. 1(f)), with a 
concomitant jump in the magnetization. 

Let us discuss the effect of thermal fluctuations. For 
this purpose, both Tc (Curie temperature) and Tn (Neel 
temperature) are estimated from MC simulations in 4x4 
clusters 17] (Fig. 1(g)). Our estimated T c is found to 
decrease with increasing Jaf, and in the E-AFM phase, 
Tn becomes almost constant. The E-phase can be disor- 
dered by thermal fluctuations faster upon heating than 
the FM phase. Provided that our estimated Tc in 2D 
system is related to Tn with the A-AFM phase of 3D lat- 
tices, Fig. 1(g) mimics well the experimental results for 
the changes of Tn by R-ion substitution in RMnOs 6]. 

Now let us turn our attention to the effect of dimen- 
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FIG. 3: (a) Phase diagram in the (x,Jaf) plane for layered 
manganites at A=0 obtained by analytic calculations, (b) 
Schematic view for spin structure of the C x Ei_ x -type phase 
at x=l/6. Hatched squares denote hole-rich C-type regions, 
(c) S(q) and (d) the DOS in the Ci/6E5/6-type phase. 



sionality. In Fig. 2(a), using bi-layer 4x4x2 clusters, the 
ground-state energy per site is depicted vs. Jaf at A=1.5 
both for numerical optimization and MF calcula- 
tions, with a good agreement between them. Figure 2(b) 
is the phase diagram on a 4x4x4 cubic lattice, obtained 
only by the MF approximation, since MC calculations 
in this case are too CPU time consuming (previous MF- 
MC comparisons suggest that this MF phase diagram is 
accurate). In the strong-coupling region, there occurs 
a chain of transitions from FM^A-AFM^C-AFM^G- 
AFM phases [19( , already obtained in 2 x 2 x 2 calculations 
|4j. The present result shows that size effects are small 
in undoped strongly- coupled manganites, which is intu- 
itively reasonable. Note that near A^1.6, a realistic value 
for manganites, the A-AFM phase is adjacent to the E- 
type state. This region could correspond to the actual 
situation observed in experiments for RMnOa: When 
the ionic radius of the R-sitc decreases, Tn of the A-AFM 
phase decreases as well, and eventually the E-AFM phase 
is stabilized for R=Ho @. In the weak-coupling region, 
the E-type phase is stable in a wide range of Jaf, as in 
the 2D calculation. 

Consider now the very interesting effect of light hole 
doping on the E-phase. Hole doping will be here stud- 
ied in the weak coupling limit, since the E-type phase is 
well understood at A=0. Figure 3(a) shows the ground- 
state phase diagram in the (x, Jaf) plane, obtained us- 
ing analytic calculations on 2D lattices at A=0. A re- 
markable feature of this phase diagram is the appear- 
ance of the novel C x Ei_ x phase, composed of long-period 
zigzag FM chains, antiferromagnetically coupled to each 
other (see Fig. 3(b) for x=l/6). As the doping frac- 
tion grows, it is expected that the previously reported 
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states at intermediate- and high-doping will eventually 
be reached HQ1 

Note that the curvature of the boundary between FM 
and C x Ei_ x phases is found to be negative for 0<x<0.5, 
indicating that phase separation occurs between those 
two phases. Here it is again stressed that the FM phase 
is metallic, while C x Ei_ x is insulating. The latter can 
be considered as a microscopic phase separated state, 
since the C- and E-type structure are mixed at the length 
scale of a lattice constant. This phase is expected to be 
made unstable easily near the phase boundary region, 
and to turn into a phase-separated state with a mixture 
of metallic FM and insulating C x Ei_ x clusters, which 
may induce CMR effects. 

In the C x Ei_ x phase, e g holes tend to localize in the 
C-type region (hatched squares in Fig. 3(b)), i.e., in the 
straight segment portion of the zigzag FM chain, and as a 
consequence charge-ordering of the stripe form is induced 
even in the A=0 limit |2fj| . This causes incommensurate 
peaks in the charge correlation at q=(27rx, 2ttx). The 
spin sector is also nontrivial (Fig. 3(c)) and the state is 
an insulator according to the DOS (Fig. 3(d)) |2lj . 

In summary, the extended phase diagrams of mangan- 



ites for x=0 and x>0 have been discussed. A novel E- 
AFM phase, stabilized at x=0 in the region of weak- and 
intermediate-couplings, is adjacent to both FM metallic 
and A-AFM states. The competition between E-AFM 
insulating and FM metallic phases suggests the possibil- 
ity of CMR effects even in undoped manganites. Several 
features of the E- to A-AFM transition, at least at the 
qualitative level, agrees with currently available experi- 
mental results. For the doped case, a microscopically in- 
homogeneous C x Ei_ x -AFM state is predicted. This state 
may contribute to the phase separation tendencies widely 
observed experimentally in Mn-oxides for 0<x<0.5. The 
discovery of these many "hidden" interesting phases in 
real manganites - undoped and doped - should be ac- 
tively pursued experimentally. This will open a new sub- 
branch of investigations in the already much interesting 
context of Mn oxides. 
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